Excited-state potential-energy surfaces of metal-adsorbed organic molecules from 
Linear Expansion A-Self-Consistent Field Density-Functional Theory (ASCF-DFT) 



o 

(N 

U 

< 

oo 
(N 



O 

C/3 



Ctf 



I 

c 
o 
o 



> 

in 
r^ 

o 
m 

> 

X 



Department Chemie, 



Reinhard J. MaureqJ and Karsten ReuteiOj 
Technische Universitat Miinchen, Lichtenbergstr. 4, D-85747 Garching, Germany 
(Dated: April 30, 2013) 



Accurate and efficient simulation of excited state properties is an important and much aspired 
cornerstone in the study of adsorbate dynamics on metal surfaces. To this end, the recently proposed 
linear expansion A Self-Consistent Field (leASCF) method by Gavnholt et al. [Phys. Rev. B 78, 
075441 (2008)] presents an efficient alternative to time consuming quasi-particle calculations. In this 
method the standard Kohn-Sham equations of Density- Functional Theory are solved with the con- 
straint of a non-equilibrium occupation in a region of Hilbert-space resembling gas-phase orbitals of 
the adsorbate. In this work we discuss the applicability of this method for the excited-state dynam- 
ics of metal-surface mounted organic adsorbates, specifically in the context of molecular switching. 
We present necessary advancements to allow for a consistent quality description of excited-state 
potential-energy surfaces (PESs), and illustrate the concept with the application to Azobenzene ad- 
sorbed on Ag( 111) and Au(lll) surfaces. We find that the explicit inclusion of substrate electronic 
states modifies the topologies of intra-molecular excited-state PESs of the molecule due to image 
charge and hybridization effects. While the molecule in gas phase shows a clear energetic separation 
of resonances that induce isomerization and backreaction, the surface-adsorbed molecule does not. 
The concomitant possibly simultaneous induction of both processes would lead to a significantly 
reduced switching efficiency of such a mechanism. 



INTRODUCTION 



Functional organic adsorbates are of special interest 
to surface nanotechnology. Their ability to selectively 
trigger dynamical changes in surface domains on molec- 
ular length scales opens many relevant applications in 
technology. A specifically interesting class of functional 
molecules are molecular switches[l , which can be re- 
versibly switched between two or more stable geometries. 
Switching of these molecules adsorbed on surfaces can 
occur via photo-excitation [2H5] or inelastic electron scat- 
tering events[B]-[S]. A typical problem to the design of 
such systems is the loss of photo-induced switching func- 
tion of gas-phase molecules upon adsorption to a surface, 
such as in the case of the prototypical conformational cis- 
trans (Z-E) switch Azobenzene (Ab) on coinage metal 
surfaces [9 . Very often this loss of function is related 
to overly strong coupling to the surface electronic de- 
grees of freedom. This coupling modifies ground-state 
energetics|10) and, even more importantly, is expected to 
change excited-state mechanisms and lifetimes. For the 
Ab case recent experiments have e.g. shown that for the 
supposedly minimized coupling of a derivative of Azoben- 
zene with bulky spacer groups, namely ieir-tert.-butyl- 
Azobenzene (TBA), the switching function can indeed be 
retrieved at Au(lll) [2111], while no switching is observed 
for the azobenzene and TBA on Ag(lll). 

A plethora of spectroscopical techniques gives access 
to the changes in electronic structure that underlie cor- 
responding excited-state dynamical processes. Still, of- 
ten enough the accessible observables do not permit to 
formulate a precise molecular mechanism, which in turn 
is a prerequisite to the understanding and subsequent 



design of corresponding systems. First-principles mod- 
elling techniques have in turn proven to be valuable tools 
for the investigation of such mechanistic details, but are 
challenged by the large system sizes and the necessity 
to simultaneously describe localized Molecular Orbitals 
(MOs) and the metallic surface band structure. Ab ini- 
tio quantum mechanical simulations, such as Density- 
Functional Theory (DFT)[TTJ H2 or post-Hartree-Fock 
approaches[13 , have a successful history as such tools 
in surface science and chemistry. The current state-of- 
the-art provides, in many cases, a reliable description 
of ground-state properties, including adsorption geome- 
tries, energetics as well as thermal barriers. When it 
comes to the description of spectroscopy and excited- 
state properties, quantum chemical approaches are the 
optimal choice for finite systems or isolating materials, 
where cluster approximations are possible. They are cur- 
rently not applicable to metallic systems though, where 
periodic boundary conditions are necessary to correctly 
describe the delocalized electronic structure. Applicable 
excited-state methods for this case include Time Depen- 
dent DFT (TD-DFT) Q3HTB], or many-body perturba- 
tion theory (MBPT) based methods [T7], the latter en- 
abling the description of both, ionic (GW) [IH] and neu- 
tral (Bethe-Salpeter equation, BSE) [15, lQlED] electronic 
excitations. In recent years, computational cost and ac- 
curacy of these approaches has improved tremendously. 
Nevertheless, current computer infrastructure and the 
remaining accuracy issues of applied density-functional 
approximations and self-energy descriptions render sys- 
tematic excited-state studies of large systems virtually 
intractable at the time. In this situation, a need ex- 
ists for highly efficient excited-state schemes that, while 
maybe not fully quantitative, allow for a qualitatively 



correct description of the major physical effects that gov- 
ern excited-state processes at surfaces. 

Already very early in the history of DFT, attempts 
to apply and/or generalize the method beyond Hohen- 
berg and Kohns rigorous proof[TT] and towards non- 
equilibrium excited-state properties have been under- 
taken. Many of them with the specific aim for a 
highly efficient description. The most rigorous and ma- 
jor extension was the Runge-Gross proof of a one-to- 
one correspondence between the time-dependent poten- 
tial and the time-dependent electron density[T3j. An- 
other line of development are functionals generalized to 
fractional occupation numbers |211 [22] , which has led to 
the standard DFT treatment for metallic systems [23], 
Utilizing a Lagrange multiplier formalism, Dederichs et 
aL[2~4] have shown how to construct constrained den- 
sity functionals 25j, constraining electrons into specific 
regions of space or spin channels. This very efficient 
method has been heavily utilized to describe electron- 
transfer processes [26H28] . but also surface reactions [29| 
[50] , Another very early approach is based on converg- 
ing the density with respect to non-equilibrium electron 
occupations that resemble excitations, so called Delta- 
Self-Consistent Field DFT (ASCF-DFT) (3JHIS3J- This 
approach, in different variations, has had a comeback 
in recent years due to its success on molecular charge- 
transfer excitations [54U4T)] . which are badly described by 
adiabadic linear-response (lr) TD-DFT using standard 
semi- local exchange-correlation (xc) kernels [HI 02] ■ Al- 
though in principle without any formal justification, this 
method has recently been put into context by a number 
of different works. Ziegler and coworkers identified a close 
connection to a constrained variational procedure [43H46] . 
which then provides a direct link to lrTD-DFT[H|. The- 
oretical works by G6rling[47] and Ayers et aL[48] in turn 
point towards a possible formal basis for an excited-state 
density functional and would, at least in the case of the 
exact xc-functional, justify a corresponding treatment. 

In the context of metal-surface adsorbed molecules 
an interesting extension to ASCF-DFT was put for- 
ward by Gavnholt et a/. [49 . This so-called linear ex- 
pansion ASCF (leASCF) scheme centers on resonance 
states that resemble gas-phase adsorbate orbitals, and 
enforces their occupation in the self-consistent density. 
This not only provides a well-defined constraint for intra- 
molecular excitations of the adsorbate, but also enables 
the description of photoemission and charge-transfer ex- 
citations. The method has already been successfully ap- 
plied to several smaller adsorbate systems [50H52"] and 
promises at least semi-quantitative results, while adding 
only little computational overhead to ground-state DFT 
calculations. In this work we strive to generalize this 
leASCF approach as an efficient means to the calcula- 
tion of excited-state potential-energy surfaces (PESs) for 
large metal-adsorbed molecules. We present necessary 
modifications to the method that allow the calculation 



of intra-adsorbate as well as substrate-mediated charge- 
transfer excitations within the same formalism. Imple- 
menting the method into the ultrasoft pseudopotential 
plane wave code CASTEP[55J, we initially compare the 
ability of this leASCF scheme to describe excited-state 
PESs in the gas phase, and establish the equivalence 
to simple ASCF for this limit of no molecule-substrate 
interaction. Thereafter we illustrate what kind of in- 
sights can be gained from this approach for the showcase 
Azobenzene adsorbed on extended Ag(lll) and Au(lll) 
surfaces. Already from the computed lowest-lying intra- 
molecular excitations along the two most important gas- 
phase isomerization pathways we can conclude on intrigu- 
ing surface-induced PES modifications with direct bear- 
ings on the isomerization mechanism. 



METHODS 

In the following section we briefly revisit the ASCF 
method and introduce the rationale behind leASCF as 
well as our modification and implementation of it. 

ASelf-Consistent-Field-DFT and leASCF-DFT 

Detailed descriptions of the ASCF-DFT approach have 
already been given numerous times in literature [321 1341 
1551 I5T1 159] . In the simplest case, excitations are mod- 
elled by reordering orbital occupation between states that 
mainly contribute to a transition. This changed popula- 
tion generates a modified density under which the Kohn- 
Sham (KS) equations [T^] are solved. Singlet excitations 
are modelled by changing populations within one spin 
channel, Triplet excitations by switching channels. Cer- 
tain care has to be taken to ensure the correct calculation 
of Singlet states. A generally used correction method is 
Zieglers sum rule[55]: E§ M = 2E S - E T . If the system 
does not show magnetization and the ground state is a 
Singlet state, Singlet excitations can also be calculated to 
a reasonable approximation without taking spin explic- 
itly into account. This has been shown for the case of gas 
phase Azobenzene [37], O2 on Al( 111)[29] and recently 
also for Iridium complexes [JD], the latter work also pro- 
viding a rationalization for the success of this approach. 

Corresponding constraints provide a reasonable de- 
scription of excited states that are well described as 
single-particle state-to-state transitions, although the 
variational adaptation of the KS states with respect to 
the excited-state density clearly does give additional flex- 
ibility. Definition of such single-electron excitation con- 
straints is a simple matter when molecular states can be 
clearly identified in character and are well separated spa- 
tially and energetically. This is almost always the case in 
minimum-energy structures of isolated organic molecules. 
More reactive geometries, i.e. transition-state structures, 



can already contain state degeneracies that hamper con- 
vergence. In such a case minimal smearing of the occu- 
pation constraint might enable calculation with only a 
small additional error in energy [57], In contrast, in the 
case of the excitation spectrum of molecules interacting 
with periodic structures, where degeneracies are ubiq- 
uitous, such a simple approach will strongly affect the 
character and the absolute energy of the excitation. In 
this situation one also has to distinguish between adsor- 
bates interacting with isolating surfaces and adsorbates 
on metals. In the first case, substrate states are mainly 
localized and generally exhibit strong hybridization with 
adsorbate states similar to interactions between two co- 
valently interacting molecules. Such hybridization can in 
principle completely modify the character and the energy 
of states, but will again generate states that are localized 
and can, in the best case, be identified in their character 
and occupied correspondingly. Therefore a simple ASCF 
approach should still capture the main part of the transi- 
tion. In the case of transition metal substrates, however, 
interactions are twofold. Following the Newns- Anderson 
mo del [53], chemisorbed molecules will show strong hy- 
bridization with d-bands, which modifies character, split- 
ting and energetic position of the frontier orbitals. Simul- 
taneously, there will also be a weak hybridization due to 
interaction with s- and p-bands. This broadens molecu- 
lar states and spreads their character over many bands 
in a small energy window. In such a case a simple ASCF 
approach, that occupies the band with the highest over- 
lap compared to a gas-phase molecular state, will miss 
significant parts of the transition and therefore strongly 
underestimate the change in density. 

Gavnholt et aZ. [49] have devised the linear expansion 
ASCF (leASCF) approach to specifically target such sys- 
tems. In their approach they do not just define con- 
straints on KS states, but on linear combinations of them. 
To illustrate this, let us shortly recapitulate the ground- 
state case for an isolated system (or for an extended sys- 
tem for each k-point separately) . There the effective one- 
particle KS equations read 



ground-state calculation this results in 
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where we define the KS potential Vks acting on the KS 
auxiliary wavefunctions and the KS eigenvalues e;. The 
density (in the following written in Dirac notation) on 
which the KS potential depends on, is constructed from 
the {IV'i}} via 
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for a finite system with N e being the number of electrons 
of the system, or in case of an extended system 
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with w k being the mathematical weight for each k-point 
and ep being the Fermi energy. In simple ASCF calcu- 
lations one instead constructs the density by replacing 
one of the states in the sum with another originally un- 
occupied virtual KS state. In leASCF, Gavnholt et al. 
propose to construct so-called resonance states from a 
linear combination of KS states instead of a single KS 
state 
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where \4%) denotes a pre-calculated reference KS state of 
the corresponding gas-phase adsorbate that ought to be 
occupied. The excited-state density then follows as 
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Equation [5] thus constructs a new KS state from un- 
occupied orbitals which resemble the chosen reference 
state and which are then used to construct the excited- 
state density. This approach can readily be used to 
model intra-molecular highest occupied MO (HOMO)- 
lowest unoccupied MO (LUMO) type excitations, where 
an equal number of electrons and holes are excited in 
the adsorbate, but also for adsorbate-substrate charge- 
transfer, where only a hole or an additional electron is 
enforced on the adsorbate states. In this case, the occu- 
pation of the remaining states has to be adjusted by low- 
ering or increasing the Fermi energy correspondingly in 
order to conserve the total electron number of the whole 
system. This approximation to an excitation can be jus- 
tified by the large ensemble of substrate electrons which 
occupy metal bands of very similar character, such that 
removing one such band from the density should induce 
only a minor error on the excitation energy. Summa- 
rizing the approach, those parts of the reference orbital, 
which are not yet included in the first N e -1 ground state 
occupied orbitals are constructed from a range of vir- 
tual KS states by projecting out components resembling 



this state and subsequently including them in the density. 
The ground-state KS procedure is therefore modified in 
the construction of the density (eq. [7| in every SCF step. 
Following this approach the kinetic energy of the system 
has to be corrected for the terms due to the newly added 
constraint orbital [15]. When breaking spin symmetry 
or including different positions in momentum space (k- 
point sampling) the procedure is followed independently 
for different spin channels or at different k-points. This 
approach is ideally suited for the description of inverse 
photo-emission and for diatomics on transition metal sur- 
faces and was shown to outperform spatially constrained 
DFT approaches as well as simple ASCF [49] . 



A fresh look on leASCF 

In the following we would like to generalize the leASCF 
approach in two aspects. First of all to allow for an ar- 
bitrary number of constraints constructed from arbitrary 
reference states without discriminating between occupied 
and unoccupied states. This provides a more consistent 
infrastructure for the description of intra-molecular as 
well as charge-transfer excitations, and might even open 
the application to systems very different from adsorbate- 
substrate complexes. Secondly, the approach should en- 
able the construction of excited-state PESs for large ad- 
sorbates in arbitrary geometries, while in the limit of 
infinite separation between adsorbate and substrate it 
should retrieve the simple ASCF result. In order to 
achieve this, certain conditions on the reference states 
\4> c ) have to be imposed. 

Modified Approach It is always possible to expand an 
arbitrary reference state \4> c ) in the complete space of KS 
states of the system under study, 
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while at the same time the remaining KS states have to 
be orthogonalized correspondingly 
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This leaves the subset of { I V 7 ^)} orthogonal to the sub- 
set of resonance states (IV'c)}; but destroys orthonor- 
mality for the complete set of KS states {\il>f) , IV'c)}- 
We therefore perform an additional orthonormalization 
on this whole set of KS states. This state transformation 
is done in every SCF step and yields a set of KS states on 
which a simple modification of the electron occupation, 
such as it is done in simple ASCF, yields an excited-state 



density as follows 
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where the only boundary condition on eq. 10 is that 
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In eq. ITT] the / e 's are given by the aspired constraint 
definition, while occupations f[ have to be adapted to 
conserve the electron number. Due to the modified con- 
struction of the KS states and the occupation reordering, 
there is no need for modification of the density construc- 
tion routine itself. The modified KS states and excited- 
state occupations also enter in the calculation of the ki- 
netic energy, which is therefore implicitly treated cor- 
rectly, again without need of further modification as was 
necessary in the original implementation of Gavnholt et 
al. in the GPAW package [55l 156] , Constraints are en- 
forced independently in different spin channels and at 
different k-space positions. 

This approach only differs from the simple ASCF ap- 
proach by the modification of the KS states, which corre- 
sponds to a unitary transformation and forces the result- 
ing KS solution to include the specified resonances. It 
naturally accounts for the hybridization-induced broad- 
ening of the adsorbate KS states at the surface and in 
all cases includes or removes the whole reference state. 
In contrast, a strong limitation of the method lies in hy- 
bridization effects that go beyond this. Due to the in- 
teraction of the sub-systems (molecule and surface) , hy- 
bridization of the system can already lead to ground-state 
occupations that are very different from the separated 
sub-system case. The correct treatment of electron tran- 
sitions then has to start from this occupation and transfer 
the corresponding amount of electrons effectively as we 
will discuss in the application to adsorbed Azobenzene 
below. 

Concerning the calculation of energy derivatives, this 
approach suffers from similar problems as the original im- 
plementation of leASCF does. The Hellmann-Feynman 
theorem does not hold due to the additional dependence 
of the non- variational coefficients of \<fr c ) on the positions 
of nuclei and the additional entropic contribution due to 
the excited-state population. A possible formulation of 
the herewith introduced |</> c )-'Pulay'-like terms still needs 
to be developed. 

Generating Suitable Reference States In the leASCF 
scheme an excited-state density is constructed that in- 
cludes a certain resonance state. All remaining states are 
variationally relaxed and therefore effectively screen the 
excitation in the self-consistent (sc) excited-state density. 
A question that remains is the selection of suitable refer- 
ence states 1 4> c ) from which to construct the resonances. 



DFT ASCF leASCF GS |eASCF EX 

£ Anon-sc sc non-sc sc non-sc sc non-sc sc 



FIG. 1. Schematic state diagram showing the frontier or- 
bitals in ground-state DFT, simple ASCF, leASCF with a 
ground-state reference and leASCF with an excited-state ref- 
erence. Projected reference states are in red. The excitation 
is visualized with an electron hole pair (filled and unfilled 
circles). Schematic orbital positions are shown for the non- 
self-consistent (non-sc) and the self-consistent case (sc). 



Such reference states could be molecular states of an ad- 
sorbate on a surface resembling an excitation, localized 
orbitals of a cluster cut-out that resemble a vacancy, or 
stemming from the very same system in a different elec- 
tronic state (as used in the case of Ab in gas phase be- 
low) . The choice depends very much on the definition of 
the sub-system and the excitation under study. The pro- 
jection restricts the resonance state itself to be an input 
quantity and it cannot change during the self-consistent 
solution of the KS equations. This stands in stark con- 
trast to simple ASCF where the non-self-consistent (non- 
sc) input orbitals from the ground-state calculation are 
used to construct the input density and are then iter- 
atively optimized to yield a self-consistent excited-state 
density(ASCF, cf. Fig. Q. 

In the work of Gavnholt et al. the choice of the refer- 
ence state fell on a virtual ground-state KS state of the 
gas-phase adsorbate (ASCF GS ). As shown schematically 
in column 3 of Fig. [T] this corresponds to calculating the 
sc excited-state solution, while forcing the constrained 
orbital into the non-sc (or ground-state) solution. This 
might be a valid approximation, if the molecular state 
of interest does not change strongly due to screening 
in the excited state. Particularly for the description of 
vertical excitation energies of equilibrium geometries or 
PESs of small adsorbates where a small number of de- 
grees of freedom defines the KS states, this might be a 
good choice. This was nicely shown for the calculation of 
excited-state PESs of small diatomics on transition metal 
surfaces [inHSU]. The single nuclear degree of freedom in 
this case reduces the chance of large variations in charac- 
ter and extent of the orbitals between ground-state and 
excited-state solutions. 

Notwithstanding, in many cases this approximation 
will fail, namely when the ground-state optimized orbital 
is not a good approximation to the final excited-state 



KS state. This is in fact the general case for the fron- 
tier orbitals of molecules with many degrees of freedom 
and/or extended 7r-systems in non-equilibrium geome- 
tries, and is especially true when applying standard semi- 
local exchange-correlation (xc) functionals. In the lat- 
ter case it is known from IrTD-DFT treatments that the 
qualitatively wrong description of ground-state molecular 
resonances in semi-local functionals hampers the descrip- 
tion of non-equilibrium geometries and charge-transfer 
excitations |37l 1571 I58j . This problem can to some extent 
be resolved by including the correct 1/r density-density- 
response behaviour into the xc-functional description[55j. 
A big strength of the simple ASCF approach is in this 
respect its additional flexibility due to the variational 
optimization of the orbitals. Although the definition of 
the excitation itself is primitive compared to IrTD-DFT, 
the additional variation allows for a consistent-quality 
description for large portions of PESs and the qualita- 
tively correct description of charge-transfer states and 
other problematic cases already with a semi-local or hy- 
brid xc-functional [5B1 I3T1 131)1 HP] . The absolute excitation 
energies will nonetheless be determined by the quality 
of the underlying xc functional, meaning that an under- 
estimation of e.g. the HOMO-LUMO gap due to self- 
interaction error will also carry over to the excited-state 
description. 

Some of this ASCF flexibility is lost due to the pro- 
jection inherent in leASCF. In order to also ensure a 
correct sc treatment of the actually constraint orbitals, 
one has to provide reference states |</> c ) that are already 
optimized to the specific excited state of interest. This 
can for example be done by calculating the simple ASCF 
solution of the excited-state reference system (here the 
gas-phase molecule) and then including reference states 
into the leASCF calculation that are already in the final 
excited state (leASCF EX , cf. Fig. [l). The solution of 
this approach would in this case correspond to the simple 
ASCF solution in the limit of zero hybridization. Such a 
generalized final-state leASCF (leASCF EX ) approach to 
arbitrary systems would thus include the following steps: 

• Calculate the electronic ground-state of the system 
of interest with DFT 

• Calculate the excited state of interest in the refer- 
ence (sub-)system using simple ASCF-DFT 

• Calculate the excited state of the system of interest 
using the excited reference state and leASCF-DFT 



Computational Details 

The method described above has been implemented in 
the ultrasoft pseudopotential plane-wave code CASTEP 
6.0.1 [S3]. The implementation for the ASCF scheme con- 
structs the changed set of KS states after every diagonal- 



ization step in the SCF procedure and uses a modified 
Fermi distribution to assign the constraint occupations, 
adapt the remaining occupations (/,') and construct the 
density from it. The newly constructed resonance KS 
state replaces the former KS state showing the highest 
overlap with the reference state. Calculations employing 
the leASCF method as implemented in CASTEP need to 
be checked for convergence with respect to the standard 
parameters of plane wave calculations such as plane wave 
cutoff and k-point sampling, but also with respect to the 
number of additional virtual orbitals that are explicitly 
included in the calculation in order to assure convergence 
of the projections from eq. [8] Standard DFT conver- 
gence enhancement methods [23, 60, 61j are used for the 
evaluation of the self-consistent density. In addition to 
the modified ASCF scheme, simple ASCF for gas-phase 
molecules in a supercell approach has been implemented. 
This is used for comparison and construction of appro- 
priate excited-state KS reference states. The implemen- 
tation of the projections in eq. [8] allowed us to use them 
also for the calculation of Molecular Orbital projected 
Density-of-States (MolPDOS) following the explanations 
of McNellis et al. [S2j. MolPDOS coefficients correspond- 
ing to gas-phase reference KS states can be printed out 
and post-processed for visualization (as shown in ref. I10[) . 
These coefficients give access to the MO occupations that 
are shown in Fig. 3 and 4. 

Excitations in this work have been modelled by effec- 
tive addition or removal of one electron in the frontier 
molecular orbitals, namely, the second highest occupied 
molecular orbital (HOMO-1), HOMO and LUMO of the 
molecule in order to describe neutral intra-molecular ex- 
citations. 

During this work we have used standard-library ul- 
trasoft pseudopotentials (USPPs)[63]. We use the xc- 
functional due to Perdew, Burke and Enzerhof (PBE) 
[6"4] throughout. Isolated Ab benchmark calculations be- 
low have been run in a 20x20x20 A supercell with a 
plane- wave energy cut-off of 350 eV and T-point sam- 
pling. The corresponding geometries along the gas-phase 
isomerization pathways have been taken from a previ- 
ous study [57], Metal-surface adsorbed Ab calculations 
were done in (111) (3x6) 4-layer surface slabs of Ag(lll) 
and Au(lll), employing a plane wave cutoff of 350 eV 
and 16 irreducible k-points. The vacuum region was 
chosen to exceed 20 A. This calculational setup closely 
follows refs [TU] and [551 where careful convergence tests 
have already been detailed. Optimized structures for the 
minimum-energy paths at the surface were taken from 
a recent work[T0]. A semi- empirical dispersion correc- 
tion was employed to ensure correct description of the 
adsorbate geometry [55J |55], although explicitly exclud- 
ing lateral interactions in order to describe the system in 
a low-coverage limit. 



RESULTS AND DISCUSSION 

In this section we apply the proposed method first to 
the isomerization of the prototypical molecular switch 
Azobenzene in gas phase and then when adsorbed on 
coinage metal surfaces. 

Isomerization of gas-phase azobenzene 

The purpose of this subsection is to demonstrate the 
equivalence of simple ASCF and leASCF in the limit 
of photo-induced E-Z-isomerization of gas-phase Azoben- 
zene. Azobenzene can efficiently Z— >E or E— >Z isomerize 
upon UV-excitation to the Sl(n ->• 7r*) and S2(7r -> ir*) 
state, respectively [57] ■ Several recent experimental [5BT- 
172] and theoretical works [751 [73H75] suggest that the 
mechanisms for both directions of isomerization mainly 
follow SI dynamics, due to strong initial population 
transfer from the S2 state to the SI. This domi- 
nant SI dynamics is very often discussed in terms of 
two mechanisms, namely, rotation around the central 
CNNC dihedral angle ui and inversion around one of 
the two NNC angles a (cf. insets in Fig. [2]), and 
there is a long-lasting controversy as to the prevalence 
of one of these two mechanisms |79l 150] . Several re- 
cent ab initio non-adiabatic dynamics studies have come 
to the common conclusion though that the dominant 
part of the photo-isomerization follows a rotation-based 
mechanism[77j ED [82] . 

In precending work we demonstrated that simple 
ASCF-DFT calculations yield a qualitatively correct de- 
scription of the involved excited-state PESs. They were 
shown to be in good agreement with higher level compu- 
tations (CC2[83]) and therefore provide a realistic repre- 
sentation of the mechanisms 3 7J. Figure [2] reproduces 
these ASCF curves (gray filled circles) following the 
ground-state optimized paths along the two main degrees 
of freedom, namely rotation around the central dihedral 
angle and inversion along one of the two central bend- 
ing angles. To summarize the picture as arising from 
these simple ASCF calculations, Z-Ab is predicted to be 
0.58 eV less stable than E-Ab, which compares nicely to 
the 0.6 eV from experiment [84] . Both pathways show sig- 
nificant barriers in the ground state, whereas only in the 
rotational pathway a low lying state-crossing with the 
first excited state can be found at mid-rotation. The 
first excited state shows a minimum at mid-rotation and 
a very small barrier at mid-inversion. The second excited 
state, in fact corresponding to a number of 7r — > n* states, 
cxibits two minima close to the positions of the ground 
state. These are on both pathways separated by sizable 
barriers. A state crossing close to the E-Ab equlibrium 
geometry with the SI state can be found. 

Figure[2]also includes the data obtained when applying 
the leASCF approach as described in the previous sec- 
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tion and using the ground-state orbitals of the isolated 
gas-phase molecule as reference orbitals at every position 
along the two pathways. Already from visual inspection 
it is possible to identify regions on both pathways where 
the difference to ASCF is minimal and regions where the 
topology is not reproduced correctly. The assumption 
that the constraint states do not change significantly due 
to the excitation seems sufficiently justified very close to 
the equilibrium geometries, but fails at the transition- 
state geometries on the S2 state, fn other words, in PES 
regions where ground-state orbitals are very good ap- 
proximations to excited-state ones the difference is min- 
imal. In contrast, in regions where due to excitation the 
orbital character and orbital ordering changes, effects can 
be quite large. In this respect, it is intriguing to note 
that the region of biggest error, namely the S2 state at 
mid-inversion, is also not correctly reproduced by lrTD- 
DFT when using ground-state PBE orbitals as a starting 
point [37] . Both effects have the same source, namely 
that GGA-DFT derived effective one-particle states are 
bad approximations to molecular resonances of the inter- 
acting many-particle system. This is especially true for 
virtual states [551155], 

Also shown in Fig. [2] are the curves calculated with 
leASCF when employing reference orbitals that were cal- 
culated with our leASCF EX approach. The correspond- 
ing results exactly reproduce the standard ASCF curves, 
because they now include relaxation effects for all KS 
states. This nicely underscores the importance of in- 
cluding state relaxation in order to generate consistent- 
quality PESs. Having established the equivalence of the 
two methods for the gas-phase limit, we now proceed in 



the next section by analyzing the effects of a metal sur- 
face on the lowest-lying excited states of Azobenzene. 



Excited-state PESs of Azobenzene adsorbed on 
Ag(lll) and Au(lll) 

Current knowledge of the mechanisms underlying 
photo- isomerization of Azobenzene and its derivatives on 
coinage metals is sparse. Photo-induced switching was 
hitherto only achieved for TBA on a Au(lll) surface, 
with a significantly reduced cross section compared to 
the gas-phase or solvent case^J[55]. The modified photo- 
absorbance of the adsorbed molecule led to the conclu- 
sion of a changed isomerization mechanism, where excita- 
tion happens indirectly via hole generation in the valence 
band and subsequent charge transfer from the molecule 
to the surface [87], At Au(lll) molecular switching was 
also achieved by resonant tunneling for Azobenzene [SUHJj 
whereas to our knowledge no switching, neither light- 
nor current-induced, has been observed for Azobenzene 
or TBA on Ag(lll). This lack of function on Ag(lll) 
surfaces strongly supports recent DFT results pointing 
to an effective loss of bistability of both derivatives on 
Ag(lll)HD|. 

Figures [3] and [4] reproduce the corresponding ground- 
state paths of Azobenzene on Ag(lll) and Au(lll) fol- 
lowing rotation and inversion. The barrier along inver- 
sion is almost unchanged, whereas the rotational bar- 
rier is strongly modified. Compared to the gas-phase 
case, the stability of the Z-Ab isomer is drastically re- 
duced from a basin depth of leV to 0.05 eV or 0.38 eV 
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FIG. 3. Upper panels: Minimum-energy paths of Azobenzene adsorbed on Ag(lll) following rotation (left) or inversion 
(right). Shown are the ground-state energy (black), the first (SI, red) and second (S2, blue) excited states as well as the 
corresponding gas-phase potential energy curves (in gray). Regions marked with dashes are of increased inaccuracy due to 
methodological restrictions further outlined in the text. Vertical dashed and dotted lines on the sides depict the position of 
E-Ab and Z-Ab minima for the adsorbed molecule. Lower panels: For both degrees of freedom, rotation and inversion, the 
integrated occupation of the projected gas-phase HOMO (dashed line) and LUMO (dotted line) in the ground state are shown. 
The horizontal line marks half filling of an orbital. 
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at Ag(lll) and Au(lll), respectively (zero-point-energy 
corrected values from ref. IT0|) . In the context of photo- 
induced E— >Z isomerization, this implies that vibra- 
tionally hot molecules on Ag(lll) after deexcitation are 
liable to thermal re-isomerization to the E-Ab isomer. 
The bottom panels of Figs. [3] and [4] show the integrated 
ground-state occupancies of the projected gas-phase ref- 
erence orbitals corresponding to the HOMO and LUMO 



of Azobenzene. For E-Ab and Z-Ab as well as follow- 
ing geometries along the inversion degree of freedom no 
considerable charge is added to or withdrawn from these 
frontier orbitals on Au(lll), cf. Fig. El on the right. 
This indicates that the bonding in these molecular ge- 
ometries is mainly physisorptive. Following the inversion 
isomerization of Azobenzene on Ag(lll), cf. Fig. [3] on 
the right, we obtain a very similar picture, although the 



Z-Ab isomer already shows some charge transfer in the 
ground state. In contrast, following rotation we see that 
on both surfaces around mid-rotation the LUMO is more 
than half occupied and the HOMO loses considerable oc- 
cupation. This is due to an orbital degeneracy of HOMO 
and LUMO at this point, which exists independent of 
metal surface adsorption. This leads to the formation of 
a strongly chemisorbed species at this point, further ra- 
tionalising the ground state barrier reduction. The shift 
of the Z-Ab minimum towards higher u angles (44°) on 
Ag(lll) together with the significant population of the 
LUMO creates a somewhat chemisorbed species in that 
case as well. 

The leASCF EX method calculates the first and sec- 
ond excited states of these surface systems by adding 
an electron to the region of Hilbert space corresponding 
to the Ab gas-phase LUMO and removing an electron 
from HOMO or HOMO-1, respectively. The correspond- 
ing excited-state curves, cf. Figs. [3] and [IJ for both de- 
grees of freedom show very similar overall topologies com- 
pared to the respective gas-phase case. When following 
inversion on Ag(lll) and Au(lll), the SI state is almost 
unchanged in comparison to gas-phase Azobenzene. A 
significant lowering of the excitation energy occurs only 
for geometries close to the Z-Ab minimum. For rotational 
isomerization, SI state energies around mid-rotation are 
reduced simultaneously with the barrier reduction in the 
ground state, while excitations close to the equilibrium 
geometries are almost unchanged. The systematic down- 
shift of the S2 state corresponds to a shift of about 1 
cV on both coinage metal surfaces all along the pathway, 
except around the mainly physisorbed E-Ab geometry. 
Two very important features for the isomerization mech- 
anism in gas phase are the state-crossings between SO 
and SI at mid-rotation and between SI and S2 close to 
the E-Ab minimum. Both can, in principle, still be ob- 
served, suggesting that an intra-molecular isomerization 
mechanism analogous to the gas phase could also prevail 
at the surface. 

As most intriguing features of surface adsorption we 
thus see a stronger lowering of the S2 state compared to 
SI, and a stronger lowering of excitations for all geome- 
tries away from the E-Ab equilibrium structure. Both ef- 
fects can be rationalized by the interaction of the molecu- 
lar dipole with the image charge that is induced in the un- 
derlying metal substrate during adsorption. Azobenzene 
in the planar trans configuration shows no significant 
dipole orthogonal to the surface in the ground and both 
excited states. Yet, following the isomerization pathways 
towards the non-planar Z-Ab isomer, the z-component of 
the dipole in the ground state increases significantly to 
a gas-phase value of 3.0 Debye (D). The corresponding 
excited-state dipole moments for gas-phase Z-Ab are 2.3 
D and 4.3 D for SI and S2, respectively. The stronger po- 
larisation of the S2 excited state thus leads to a stronger 
interaction of the molecular dipole with the image charge 



TABLE I. Vertical excitation energies (in eV) for the first 
Singlet excited state of Z-Ab (SI) and the second Singlet ex- 
cited state of E-Ab (S2), as well as the difference between 
them. Reference values shown are taken from experiments in 
solvent |67l 188] and from high-level quantum chemical (RI- 
CC2) calculations for the isolated molecule |37| . 





Z-Ab SI 


E-Ab S2 


Energy Difference 


Exp. a 


2.87/2.92 


3.89/4.12 


1.02/1.2 


CC2 ©gas 6 


3.00 


4.07 


1.07 


ASCF-B3LYP @gas 6 


2.30 


3.33 


1.03 


ASCF-PBE @gas 6 


2.10 


2.98 


0.88 


leASCF-PBE @Ag(lll) 


2.03 


2.44 


0.41 


leASCF-PBE @Au(lll) 


2.27 


2.38 


0.11 



refs. [67] and [ 



ref. [37] 



monopole and explains the particularly pronounced low- 
ering of the S2 PES upon adsorption obtained in the 
lcASCF EX calculations. An important point to mention 
here is that the variational treatment in leASCF and 
ASCF approaches enables such an image charge build-up 
due to polarisation effects (opposed to IrTD-DFT treat- 
ments), although we emphasize that this is unlikely a 
quantitative account. 

Another effect that modifies excited-state behaviour 
is the hybridization of molecular with surface states. A 
marker for the strength of hybridization is the change in 
ground-state occupation of the frontier orbitals, which we 
find much more pronounced for Ab adsorbed on Ag(lll) 
than on Au(lll). In regions where orbitals show oc- 
cupancies very different from the gas phase, e.g. Z- 
Ab at Ag(lll), we obtain PES changes that are more 
significant than for regions where occupancies do not 
change drastically. This effect is especially strong around 
mid-rotation, where the ground-state occupation of the 
LUMO already increases beyond one electron on both 
surfaces. This prohibits the full transfer of one further 
electron into the LUMO in the leASCF EX excited-state 
calculations and we instead only perform these calcula- 
tions by enforcing a full two-electron occupation of the 
LUMO. In Figs. [3] and [4] we mark these regions with 
dashed lines to emphasize the expected increased uncer- 
tainty due to the concomitant violation of the excitation 
constraint. We believe that these parts of the SI curves 
can only serve as an upper estimate to the actual PES 
topology and attest that such situations of strong hy- 
bridization and charge transfer represent a clear limita- 
tion to the leASCF approach. 

Notwithstanding, even when only taking them qual- 
itatively, the obtained results clearly show that an ex- 
plicit treatment of hybridization, charge transfer and im- 
age charge effects is necessary to appropriately describe 
ground- and excited-state PESs of a functional molecule 
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like Azobenzene when adsorbed at metal surfaces. In- 
vestigating Ab in a van-der-Waals potential to merely 
mimick the effect of surface-modified molecular geome- 
tries on the switching function, Flofi et al. [5S] recently 
reported only a small increase in conversion times and 
decrease of photo-yield compared to the gas phase, while 
otherwise the photo-isomerization was unaffected. With- 
out yet embarking on actual dynamical simulations, the 
leASCF-obtained PESs topologies already indicate that 
much larger effects are induced by the metal electronic 
structure. Notably this is the image charge induced low- 
ering of excited-state PESs, which due to the varying de- 
gree of state polarization and molecular dipole moment 
does not occur globally, but differentially "skews" indi- 
vidual state topologies and vertical excitation energies. 
For the present system this leads to a strong lowering of 
the S2 state particularly around the Z-Ab geometry. 

By itself this image-charge lowering might already ra- 
tionalize a significantly reduced switching efficiency of 
the traditional intra-molecular gas-phase isomerization 
mechanism at the surface: This mechanism proceeds via 
initial excitation to S2 and fast population transfer to SI 
for the E— >-Z isomerization, as transition to SI is symme- 
try forbidden in the E-Ab geometry. The back-reaction 
Z— s-E instead involves direct excitation to SI. In the gas 
phase the corresponding vertical excitation energies for 
the two reactions differ substantially (cf. Table |l| and 
allow the two isomerizations to be selectively induced by 
light with two largely differing wave lengths. In contrast, 
at the surface our leASCF EX results suggest that the 
selective image-charge induced S2 lowering reduces this 
difference for the two transitions substantially. While in 
the gas phase it amounts to more than 1 eV, particularly 
at Au(lll) the difference between E-Ab S2 and Z-Ab 
SI reduces to 0.11 eV. Considering an additional state 
broadening at the surface, this proximity of the two dif- 
ferent excitations alone might then already cause a signif- 
icant loss of switching efficiency via this intra-molecular 
mechanism as the forward and backward isomerization 
can simply no longer be selectively triggered. 

At least qualitatively, these findings should also be ro- 
bust against the other clear limitation of leASCF, namely 
the one imposed by the employed approximate DFT func- 
tional. Already in the gas phase our preceding work 
demonstrated that GGA-PBE based ASCF (but also 
GGA-PBE based IrTD-DFT) severely underestimated 
absolute vertical excitation energies for Azobenzene com- 
pared to accurate quantum-chemical (PJ-CC2) calcula- 
tions [37J . These were primarily global shifts of entire re- 
spective excited-state PESs though and largely left topo- 
logical features like barriers unchanged. Addition of ex- 
act exchange can remedy these self-interaction induced 
shortcomings of the semi-local functional for gas-phase 
Azobenzene [37]. However, simultaneously it would re- 
move much of the balanced error cancellation in the 
description of the metal substrate |90[ |9"T] . For metal- 



surface adsorption there is at present no feasible and 
equally efficient alternative to semi-local DFT. GGA- 
PBE based leASCF excited-state PESs thus have to be 
seen in light of the self-interaction induced overpolariz- 
ability and wrong relative positions of molecular and sub- 
strate states, which will affect the observed image charge 
and hybridization effects upon adsorption. While thus 
certainly not quantitative, the approach still enables in 
our view an effective first account of the electronic struc- 
ture and charge distributions at the metal surface and 
is thus a viable means to generate further insight into 
the intricacies of surface functionality of large organic 
molecules like Azobenzene. 



CONCLUSIONS 

We presented an alternative implementation of the 
leASCF method of Gavnholt et al. [49] and neces- 
sary modifications to allow its application to complex 
metal-surface adsorbed chemical reactions. The current 
method provides a computationally efficient way to de- 
scribe low-lying localized excited states in large periodic 
systems. The correct calculation of reference states that 
are used to generate the resonances assures consistent 
quality of excited-state potential energy surfaces and also 
sets the connection to simple ASCF in the limit of van- 
ishing hybridization, in the surface context between ad- 
sorbate and metal substrate. For the example of photo- 
induced isomerization of Azobenzene at (111) coinage 
metal surfaces we illustrated that the approach yields an 
account of the additional stabilization due the interac- 
tion of large excited-state dipoles with the substrate im- 
age charge and hybridization-induced state renormaliza- 
tion. As such the method at least qualitatively describes 
most important physical effects that arise from the in- 
teraction with the electronic structure and charge distri- 
butions at the metal surface, and thus allows to discuss 
surface effects on the molecular functionality beyond the 
level of surface-modified adsorbate geometries. For the 
Azobenzene showcase system this is already highlighted 
by the observed unbalanced shifts of the intra-molecular 
excited-state energies of E-Ab and Z-Ab, which leads to 
an alignment of isomerization-inducing resonances. Re- 
ducing the ability to selectively trigger back and forth 
isomerization, this could be a first important piece in the 
puzzle to understand the strongly reduced isomerization 
efficiency at the surface - at least of the traditional gas- 
phase mechanism. 

We believe that the approach presented in this 
work, although approximate in nature, enables a semi- 
quantitative account of excited-state properties for large- 
scale systems and might prove to be very useful specifi- 
cally for large hybrid organic/metallic interfaces. While 
it may never replace theoretically rigorous methods, such 
as TD-DFT or many-body perturbation theory, it fills a 
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gap in the current methodological spectrum, where these 
more accurate methods are not yet applicable due to their 
computational expense or where currently used approx- 
imations in the xc-kernel or self-energy description in 
these methods cause a lack of consistent accuracy. In- 
dependent from the development of these schemes there 
will always be the need for very efficient treatments that 
allow fast screening on a qualitative or, when solid bench- 
marking is done, possibly semi-quantitative level. 
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